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ABSTRACT 





The present paper an attempt to study flow of an unsteady dusty fluid through a circular pipe under the influence of impulsive pressure gradient. We 
have assumed that the initial fluid velocity to be some average velocity. Analytical expressions for the velocities of fluid and dust particles are 
obtained by using Laplace transform technique. The influence of mass concentration and relaxation time of the dust on the velocity have been 
clearly pointed out. Results have been discussed with the help of graphs. 
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1. INTRODUCTION: 

The mechanical behaviour of dusty fluids has been a subject of study 
receiving greater attention during the recent past of several researchers 
in the field of Fluid Dynamics. Problems dealing with the influence of 
dust particles on viscous flows find place in several branches of 
Science and Technology. Some such flows are those of red corpuscles 
and other bodies in blood, fluidization, combustion, environmental 
pollution and chemical engineering flow in rocket tubes where small 
carbon or metallic fuel particles are present and so on. 


Since a dusty fluid is a mixture of fluid and fine dust particles, its study 
is of practical importance. Particularly, the study of flow of dusty 
fluids through circular pipes is of great importance which has 
applications in various Industries. Considering blood as a two-phase 
fluid, this study also gives some insight into the blood flow through 
veins or arteries 


Interest in problems of mechanics of systems with more than one 
phase has developed rapidly in recent years. Much work has already 
been done on dusty fluid flow. The equation n of motion of fluid 
carrying small particles were first derived. P.G.Saffman (1962) who 
investigated the stability of the laminar flow of a suspension in which 
the particles are uniformly distributed. A.H. Nayfeh (1966) introduced 
equation of motion of fluid containing dust particles taking volume 
fraction of dust particles into account. S.Uchida (1956) has made some 
theoretical analysis of pulsating viscous flow in circular pipe. S.K.Nag 
et al., (1981) have considered the flow through an elastic tube rather 
than a rigid tube. Kishore and Pandey (1977) have discussed the same 
problem but they found the effect of particles size on velocity of 
sedimentation. Datta and Dalal (1992) have discussed the unsteady 
flow of dusty fluid through a circular pipe. S.N. Dube and Srivastava 
(1972) have discussed the unsteady flow of dusty viscous fluid in a 
channel and pipe under the influence of linear pressure gradient. 
Sambasiva Rao (1969) has obtained the analytical solutions for the 
dusty fluid flow through a circular tube under the influence of constant 
pressure gradient. S. Rashmi et.al. (2007) have discussed the unsteady 
flow of a dusty fluid between two oscillating plates under varying 
constant pressure gradient and Siddapasappa et.al. (2008) have 
discussed the viscous dusty fluid flow with constant velocity 
magnitude. 


In all the above discussed papers the authors considered the zero initial 
velocity, but in this paper, the authors study the problem of unsteady 
flow of dusty fluid through a circular pipe induced by an impulsive 
pressure gradient taking into account the volume fraction of the dust 


particles and we have assumed that the initial fluid velocity to be some 
average velocity. Analytical expressions have been obtained for the 
velocity components of both of the dust and the fluid. The influence 
of mass concentration and relaxation time of the dust on the velocity 
have been clearly pointed out. This problem finds its application in the 
blood flow through Arteries. 


2. OBJECTIVES: 
1. To know the behavior of blood flow 


2. To know the effect of time on velocity of fluid 


3. METHODOLOGY: 
The Laplace transform method is use in this research paper and result 
discuss through graphs using Origion software. 


4. BASIC EQUATIONS OF MOTION AND THEIR SOLUTION: 
Consider the flow of in incompressible dusty fluid through a circular 
pipe of radius a. Let the flow be induced by an impulsive pressure 
gradient along the axis of the pipe taken as x — axis. The equations of 
motion of unsteady viscous fluid with uniform distribution of dust 
particles are given by (Saffman, 1962) 
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Where U ,V are the components of the fluid and particle velocity 


in the direction of x- axis which is taken along the length of the tube. 
The last term on the right hand side of equation (1) represents the force 
exerted by the particles on the flow while the term on the right hand 
side of equation (2) is a similar force term exerted by the fluid on the 


particles. {= Nm is the ratio of the mass density of the particles 
p 

and the fluid, commonly known as mass concentration of the dust 

particles. G = Sad t= I is the relaxation time of the particles. 


is the fluid density, / is the viscosity of fluid, N is constant number 
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density of the particles, @ is volume fraction of particle phase, D 


is the pressure of fluid. K is stokes resistance coefficient. (i.e. K = 3 


*k 
Li 7 d, dis particle diameter) and m is mass of each particle, 7 is 
the radial distance. 


Initial and boundary conditions are 
=U,v =U at <0,0<r <a (3) 


u=Oat >0, r =a (4) 


Introducing the following dimensionless variables: 
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Where U being a reference velocity and Py constant hydrostatic 


pressure, Equations (1) and (2) reduce to 
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Where F(t)= i op , R= 2aUp _ Reynolds Number, 
2 ax lu 
f= Nm The mass concentration of dust particle (ratio of the mass 
fe 


density of the particles and the fluid), B= ae ( —— - particle 
Ur k 


The pertinent initial and boundary conditions are 


relaxation time 


u=v=lat<00<r<l (7) 


u=Oat>O,r=l (8) 


Applying Laplace transform with respect to time t and using the initial 
conditions given in Equation (7), Equations (5) and (6) reduce to 
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Where {@(s), ¥(s), F(s)}= [{u(t), v(t), Fe dt 
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From (10), we get 
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v(s) = u + (11) 
s+ Be s+ B 
Substituting the value of Vin (9), we get 
2 ~ = 
AL eS (12) 
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Solving Eq. (12) satisfying boundary condition U = Oat r = 1 and 
the regularity condition at r = 0 we get 


[R 
x oo gar ag 
FF 2 


Tie ee (13) 
bb R 
PARA 
2 


Using the convolution theorem and Cauchy’s residue theorem, inverse 
transform of Eq. (13) can be written as 
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A. ’s are the zeros of the Bessel function sf 2) i.e 


b R= -22, * hence b _(s, + B)= fB s,+s,(1- 
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= oT 


Applying inverse Laplace transform to Eq. (11), we get 
v(t) = exp(—B t) + Bl u(r) exp(— Bit- t)) dt (16) 
0 


5. IMPULSIVE PRESSURE GRADIENT CASE: 

Consider a particular case of the impulsive pressure gradient in the 
form F(t) = P H (t ), where P represents the constant adverse pressure 
gradient and H (t) is Heaviside function of time t. 


For this case, the fluid velocity given in Eq. (14) can be written as 


u(t) = ie —1)+ 


(17) 
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Substituting (17) in (16), we get 





v(t) = exp(— “B+ ade it) —exp(-ft) (1—exp(— m), BAD, [seeder A) 


s, +B B S, 5, +B 
(18) 
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6. INSTANTANEOUS RATE OF DISCHARGE: 
The instantaneous rate of discharge of fluid is given by 


* 1 
Q,= g 2| urdr (19) 
0 


S 


Similarly the dimensionless instantaneous rate of discharge for 
particle is 


1 

QO, = 2| vrdr 
0 

(20) 


Using the values of u and v from Equations (17) and (18) in Eqs. (19) 
and (20) respectively, the time dependent rate of discharge of fluid and 
particle phase are obtained. 


7. RESULT AND DISCUSSION: 

In order to discuss results, numerical computation for Equations (17) 
and (18) has been made on using typical values of the dimensionless 
parameters e.g. volume fraction 0, mass concentration f, adverse 
pressure gradient P, Reynolds number R, and for various values of 
time t. In order to compute the series solution for u- fluid velocity and 
v — particle velocity, we have considered the first 275 zeros of Bessel 
function of Zeroth order and first kind numerically. Considering the 


density ration Ps O(10) and the parameter B = 1.0, it is noted that 


f 


g@ = —. The fluid velocity and rate of discharge at different times 
10 

have been computed for various values of the parameter f, (i.e. @ also) 

and they are presented in Figures. la — 4a. 


Figurela shows the velocity profiles of fluid and particle phase for 
various values of f. From the graphs we find that the velocity of the 
fluid as well as that of the particle phase decreases with the increase 
of mass concentration f. i.e. of particle loading and volume fraction 9. 
Comparison of the graphs shown in Fig.la revealed that the velocity 
of the particle phase is always less than that of the fluid phase, but near 
the wall velocities of particle and the fluid phase approach to be equal 
to each other as f increases. It 1s also noted that velocities of both fluid 
and dust particles are maximum on the axis of the tube and fluid 
particles moves faster than the dust particles everywhere in the tube. 


Figure 2a shows the distribution of velocity of the fluid and that of the 
particle phase at the center of the pipe with respect to f for different 
values of t. The difference between the velocities of fluid and particle 
phase tends to zero as time increases and they are nearly equal at for 


t=100. 


Figure 3a gives the graphs of the velocity of fluid and particle phase 
against time t for various values of particle loading f. From figure it is 
observed that the velocities of fluid and particle phase first increases 
with time due to acceleration up to a constant value for large time as f 
increases and the steady state 1s reached at ttme ¢ >100 for both fluid 
and particle phase. 


Figure 4a shows the growth of discharge rates of fluid and particle 
phase with time. It is seen that these rates of discharge decreases first 
with time and increase with increase of f (or ). For large time i.e. 
t >100 the discharge rates attain the steady state due to the steadiness 
of velocity for both fluid and particle phase. Further, the time to reach 
the steady state increases with the increase of f. 


Also we compare these results with the results obtained by Datta and 
Dalal [10] and they are presented in Figures. | — 4. We observed the 


changes in fluid and particle phases and also the discharge rate of fluid 
and particle phase between 20 < t <100. 
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Fig4a. Variation of discharge rate of fluid and particle phase 
with time when p=0.1, R= 150 6=1.0 
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